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ABSTRACT 


The  activites  at  the  ETS  have  expanded  from  observations  of 
deep-space  artificial  satellites  during  the  dark  hours  to 
observations  of  near-Earth  satellites  near  noon.  Not  only  do 
these  varying  tasks  require  different  video  cameras  and 
associated  equipment,  they  can  benefit  from  a  tailored  analysis 
of  the  resulting  astrometric  data.  In  particular,  the 
observation  of  near-Earth  artificial  satellites  very  near  the 
astronomical  zenith  has  become  a  frequent  mode  of  operation.  For 
some  of  these  satellites  a  parallax  can  be  discerned  when  both 
ETS  telescopes  are  used.  For  others,  their  nearly  circular  orbit 
coupled  with  a  zenith  passage  suffices  to  fix  their  element  set. 
This  Project  Report  details  both  types  of  analysis. 
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I. 


INTRODUCTION 


Certain  artificial  satellites  are  best  illuminated  at 
twilight.  The  reason  is  they  are  so  close  to  the  Earth  that  were 
it  darker  they  would  be  eclipsed.  Clearly  were  it  brighter  than 
twilight,  their  phase  angle  would  be  larger  and  the  sky 
background  would  be  very  bright.  A  consequence  of  these 
considerations  is  that  they  are  most  easily  observed  near  the 
zenith.  A  result  of  their  closeness  to  the  Earth  is  increased 
atmospheric  drag.  Therefore,  whether  by  design  or  owing  to 
natural  forces,  their  orbits  are  nearly  circular.  Hence,  there 
are  only  four  meaningful  entries  in  their  orbital  element  set, 
not  six. 

The  bulk  of  this  Report  is  concerned  with  two  different 
methods  of  obtaining  these  four  quantities  as  a  result  of 
observing  the  nearly  zenithal  passage  of  an  artificial  satellite 
in  a  nearly  circular  orbit.  Two  cases  are  presented.  In  the 
first  instance,  a  satellite  beyond  about  1000  km  will  move  slowly 
enough  across  the  small  field  of  view  of  the  ETS  telescopes  to 
allow  for  a  good  deduction  of  its  angular  velocity.  This  plus 
the  position  and  the  relevant  celestial  mechanics  uniquely 
determines  the  element  set.  In  the  second  case  the  satellite  is 
much  closer  and  moves  too  fast  to  reliably  determine  its  speed. 
However,  its  very  closeness  means  that  it  has  a  measurable 
parallax.  The  measurement  of  the  parallax,  the  slope  of  the 


streaks,  and  the  position  also  uniquely  defines  the  element  set. 

This  Project  Report  discusses  both  of  these  cases  in  depth. 
The  principal  aspect  of  the  analysis  missing  is  an  estimation  of 
the  errors  involved  for  we  usually  observe  an  object  slightly  off 
the  zenith,  whose  orbit  is  not  truly  circular,  and  so  on.  As 
observing  needs  demand,  more  sophistication  will  be  introduced. 


II.  POSITION  PLUS  ANGULAR  VELOCITY 


A.  Data  Acquisition 

From  a  videotape  of  the  passage  of  a  near-Earth  satellite 
played  back  in  very  slow  motion  one  can  clearly  see  the 
progression  of  the  head  of  the  streak.  If  the  recorder  has  a 
stop  frame  mode,  then  individual  frames  (and  fields;  the  1/30  of 
a  second  video  frame  consists  of  two  1/60  of  a  second  interleaved 
video  fields)  can  be  viewed.  The  addition  of  a  device  for 
measuring  location  on  the  videotape  and  a  frame  counter  complete 
the  data  acquisition  hardware  (see  Figs.  1  and  2). 

With  such  a  setup  (Fig.  2)  one  can  make  a  progression  of 
pairs  of  (say)  rectangular  coordinate  measures  of  the  head  of  the 
streak.  A  simple  method  of  treating  the  data  that  smooths  out 
non-linearities  in  the  camera  target,  stretching  of  the 
videotape,  and  so  on,  is  to  difference  them  pairwise.  Division 
by  the  appropriate  number  of  frames,  multiplication  by  30,  and 
then  multiplication  by  the  plate  scale  gives  the  velocity  in  arc 
seconds  per  second  of  time.  Averaging  a  series  of  such  pairwise- 
differenced  velocities  across  the  field  of  view  provides  the 
angular  velocity  data  for  analysis. 

It  is  useful  to  be  able  to  determine  the  orbital  inclination 
from  the  slope  of  the  streak.  To  this  end  the  pairs  of  points 
are  fit,  via  the  method  of  least  squares,  to  a  straight  line. 

The  last  piece  of  information  needed  is  the  position  and  we  know 
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Figure  1.  Schematic  representation  of  a  succession  of  overlaid 
video  frames  showing  the  progressive  movement  of  the  head  of  a 
satellite's  streak. 


that;  it  is  at  the  zenith  at  the  time  of  observation. 

B.  Celestial  Mechanics 

From  the  solution  of  the  two-body  problem  we  can  express  the 
geocentric  right  ascension  and  declination  in  terms  of  the 
orbital  element  set.  Using  the  symbols  defined  below,  the 
results  are 

a  -  ft  +  tan" 1 ( cositanu)  ,  6  »  sin~ 1 ( sinisinu)  (1) 

a  =  semi-major  axis 
e  -  eccentricity 
a)  *  argument  of  perigee 
i  =>  inclination 

n  =  longitude  of  the  ascending  node 

T  »  time  of  perigee  passage 

v  *  true  anomaly 

u  =  v  +  u)  =  argument  of  latitude 

E  =  eccentric  anomaly 

M  =  mean  anomaly 

n  *  mean  motion 

a.  =  geocentric  right  ascension 

6  -  geocentric  declination 

t  »  local  sidereal  time 


geocentric  latitude  of  the  observer 


A  -  topocentric  right  ascension 
A  -  topocentric  declination 

Because  I  have  assumed  the  orbit  is  circular,  e  •  0.  I  can  now 
define  T  and  Q  such  that  ai  -  0 .  Then,  at  time  t, 

u-v+u-v-M-E-  n(t  -  T) 

Rewrite  the  expression  for  the  mean  anomaly  as 

M  «  Mz  +  n(t  -  tz) 

where  tz  is  the  Universal  Time  corresponding  to  zenith  passage. 

As  I  will  explicitly  show  below,  at  the  zenith, 

a  »  A  -  t  and  6  -  A  »  <J> ' 

for  any  observer.  Thus,  since  a  =*  x  and  6  =  <j> '  at  the  zenith  by 

definition,  it  follows  that  Eqs.  (1)  implies  Eqs.  (2); 

~z  -  8  +  tan-1  (cositanMz)  ,  <p ’  -  sin" 1  ( sinisinMz)  (2) 


For  our  observatory  4>  ’  >  0.  Hence,  as  i  e  [  0 ,  tt  ]  and  sini  >  0,  it 
follows  that  Mz  e  [0,ir]  too.  Equations  (2)  are  two  equations 


in  three  unknowns.  The  angular  velocity  provides  two  more 
equations  and  introduces  the  fourth  unknown  (i.e.,  a).  Before 
seeing  this  I  need  to  project  the  geocentric  motion  onto  the 
topocentric  celestial  sphere. 

C.  Topocentric  Appearance 

Let  the  observer's  geocentric  location  be  given  by  £  - 
p_fc(-r,<J>')  where  l_  is  a  unit  vector  of  direction  cosines  in  the 
equatorial  coordinate  system.  The  satellite's  geocentric 
location  in  this  reference  frame  can  be  written  as  £  -  rSL(a,6) 
The  topocentric  location  is  R  -  R£(A,A), 

R  -  £  -  £ 

Equation  (3)  can  be  rewritten  in  component  form  as 


psinh 

tan(o  -  A)  - 


qsin(Y  -  5) 

tanU  '  A)  “  T'-  q'cTsT~r~5y 


R/r  =  sin(6  -  Y)csc(A  -  Y) 


where  the  auxilary  quantities  p,  q,  and  Y  are  defined  by 


p  -  (  p/r )  coS(|>  •  see6 


q  -  (p/r)sin4>'cscY 


(5 


tanY  -  tan<J>'cos[  (a  -  A)/2]sec[h  +  (a  -  A)/2] 


h  is  the  hour  angle  -  t  -  a. 

Note  that  h  -  0  when  t  -  a,  that  is  for  any  object  on  the 
celestial  meridian  there  is  no  parallax  in  right  ascension. 
Therefore,  a  »  A  at  the  zenith.  Note  too  that  if  h  -  0  and  a  - 
then  Y  -  <j> '  .  Whence  if  <5  -  too  (i.e.,  at  the  astronomical 

zenith),  6  -  A  from  the  middle  of  Eqs.  (4).  This  proves  my 
assertions  earlier. 

Differentiate  Eq.  (3)  once  with  respect  to  the  time  t  and 
solve  for  R,  A,  and  A.  Evaluate  the  results  at  the  zenith 
[wherein  R  -  r  -  p  despite  the  apparent  indeterminateness  in  the 
last  of  Eqs.  (4)].  You  will  find 


"  0  * 


I  have  replaced  r  by  a  because  e  is  zero. 


(6 


From  the  differentiated  version  of  Eqs.  (1)  and  the  fact 


that  u  -  n(t  -  T)  in  this  instance  one  finds 

.  ncosisec2M 

a  *  T~~co's~*Ttan*M‘  6  “  naini3ec<> •  cosM 

Now  it  is  known  that  Mz  e  [ 0 , tt  ] .  So,  either  Mz  is  an  acute 
angle  and  sinMz  -  sin<|>'csci  or  Mz  is  an  obtuse  angle  and  can 
be  written  as  Mz  -  n  -  mz  where  m2  e  [o,n/2].  As  in  this 
instance,  sinMz  -  sin(ir  -  mz)  -  sinmz  -  cscisin«j>',  tanMz 
-  ♦_  sin<|> '  /  ( sin2  i  -  sin2^')1/2  in  the  two  cases.  Only  the  square 
of  this  is  needed, 

az  -  ncosisec2  (j) '  ,  6Z  «  nsinisec<j> '  cosMz 

Observe  that  if  Mz  is  less  than  or  equal  to  tt/2  then 
cosMz  -  +(sin2i  -  si n24»')l/z  csci  whereas  if  Mz  is  greater 
than  or  equal  to  tr/2,  cosMz  -  -cosmz  -  -(sin2i  -  sin2^')1/2 
csci.  Therefore, 

•  i/2  M  acute 

<5z  -  +_  n(sin2i  -  sin2<f>')  secd>  *  z 

M  obtuse 
z 

From  Eqs.  (6)  it  follows  that  the  sign  of  Az  is  the  same  as  that 
as  6Z.  Moreover,  sgn(Az)  is  a  quantity  that  is  determined  by 
the  observations  themselves.  Whence,  the  correct  quadrant  for 
Mz  can  be  ascertained  by  observing  the  north-south  direction  of 
motion.  Finally,  for  the  topocentric  angular  velocity, 


nacosisec2<|> '  -  pt  .  +  na(sin2i-sin24>' ) 


a2  -  - 


Before  proceeding  to  a  calculation  of  the  motion  as  seen 
projected  onto  the  telescope's  focal  plane,  which  is  what  is 
recorded  by  the  videotape  recorder,  consider  the  topocentric 
angular  speed  9Z, 


8Z  -  (Azcos2Az  +  Az)1/ 


na  .  .  »/* 

- - -  [l  -  2(pi/na)eosi  +  (  pi/na)  2cos2<J> '  ] 

a  -  p 


For  a  polar  orbit,  and  many  near-Earth  satellites  are  nearly  in 
such  an  orbit, 


8z  “  a  'p'  t1  +  (  P  i/na)  2cos2()>  ’ 


na  1 

=  ~  i  ~  [l  +  2^Pi'/na)  2c°s24>’  ] 


as  na  >>  pi.  Thus,  the  binomial  theorem  expansion  is  especially 
accurate  in  this  case.  It  should  also  be  clear  that  for  a  fixed 
value  of  a,  9Z  is  a  maximum  for  a  polar  orbit.  Table  1 
provides  a  matrix  of  9Z  values  (in  deg/sec)  for  Rz  „  100(100) 
1500  km  and  i  -  35(10)85,  90°. 


D.  Focal  Plane  Coordinates 


The  motion  we  are  interested  in  appears  to  occur  on  the 
surface  of  a  sphere.  We  analyze  it  after  projection  onto  the 
focal  plane  of  the  telescope.  The  nature  of  this  projection  is 
complicated  but  for  small  fields  of  view  the  results  are  simple. 
So,  incorporating  this  new  approximation  into  our  analysis,  the 
rectangular  coordinates  5,  n  in  the  camera's  focal  plane 
(measured  in  the  same  units  as  the  telescope's  focal  length  f), 
where  n  is  positive  northward  (increasing  declination)  and  £  is 
positive  eastward  (increasing  right  ascension),  are  given  by 

5/f  -  (A  -  A* ) cosA* ,  n/f  -  A  -  A*  1*) 

In  Eqs.  (9)  5  and  n  are  the  rectangular  coordinates  appropriate 
for  an  object  near  the  optical  axis  of  the  telescope  when  the 
object's  topocentric  right  ascension  and  declination  are  A,  A. 

The  topocentric  equatorial  coordinates  of  the  point  where  the 

(extended)  optical  axis  of  the  telescope  would  pierce  the  sky  are 

.  *  .  » 

A  ,  A  . 

Now,  since  I  have  assumed  that  the  field  of  view  is  small, 
it  follows  that  for  the  cases  of  interest  (n  large)  the  transit 
time  across  the  field  of  view  will  be  short.  Therefore,  a 
sufficiently  accurate  representation  of  A  and  A  around  the  time 
of  zenith  passage  is  just  the  first  two  terms  in  a  power  series 


Az  and  Az  are  just  equal  to  t  and  <(> ' .  Az  and  Az  were 

last  seen  in  Eqs.  (7).  Combining  this  information  with  Eqs.  (9) 

and  the  fact  that  A*  -  t,  A*  -  yields 

£/f  -  Az(t  -  tz)cos(J>',  n/f  -  Az(t  -  tz) 

The  speed  of  motion  across  the  field  of  view  is  (£2  +  n2)1/2 
which  is,  it  turns  out,  just  f 9 z ;  0Z  was  given  in  Eq.  (8). 

The  slope  of  the  streak,  dn/d£,  is 

dn  na(sin*i  -  sin2$')  1  sec2$’ 

3  “  dT  ”  -  (nacosisec*^1  r  "p  t  )'  ^ 1 C 

where  the  plus  or  minus  sign  ambiguity  still  stems  from  the 
quadrant  ambiguity  for  Mz.  Take  the  limit  as  na/pi— Then, 

0Z  —  [GME/a(a  -  p ) * ] * / 2  (11 


since  n2a3  -  GMg,  and 


s  — (sin2i  -  sin24>’)1/2  seci 


Or,  solving  for  the  inclination  (s  is  the  slope  of  the  streak 
determined  by  the  least  squares  fit  referred  to  earlier) 


(1  ♦  s2)  ' 


Now  the  plus/minus  sign  ambiguity  is  from  a  square  root 
operation. 


In  this  approximation  one  can  proceed  as  follows: 

1  .  Note  the  direction  of  motion.  This  fixes  the  sign  of 
« 

Az  and  the  quadrant  for  Mz. 

2.  Calculate  f0z  from  the  videotapes  as  discussed  in 
subsection  A.  Solve  the  cubic  equation  for  the 
semi-major  axis  a  in  Eq.  (11a). 

3.  Determine  the  slope  of  the  streak  as  discussed  in 
subsection  A  and  use  it  in  Eq.  (11b).  This  provides  a 
value  for  the  inclination  i.  The  correct  quadrant  for 
i  is  fixed  by  the  second  part  of  Eqs.  (2).  Obtain  a 
value  for  Mz  from  this  expression  too. 

4.  Use  the  time  of  zenith  passage  plus  the  first  part  of 
Eqs.  (2)  to  calculate  the  longitude  of  the  ascending 
node  fi. 


Refine  a  and  i  as  necessary  from  Eqs.  (8,10)  in  an 
iterative  fashion  and  repeat  steps  3t  4,  5. 


5. 


E.  The  General  Problem 


Equations  (1)  and  their  derivatives  are  four  equations  in 
the  four  unknowns  a,  i,  £2,  and  T.  Projecting  them  onto  the 
topocentric  celestial  sphere  via  Eqs.  (4,  5)  introduces  no  new 
variables.  Actually  there  is  a  better  representation  in  terms  of 
topocentric  variables  given  by 

sin(a  -  A)  «  psinH 

sin(<5  -  A)  -  qsin(Y  -  A) 

where  p,  q,  and  Y  are  defined  in  Eqs.  (5)  and  H  -  t  -  A. 

Finally,  the  modelling  of  the  motion  across  the  focal  plane  in 
Eqs.  (9)  also  introduces  no  new  variables.  Thus,  £,  n  and  their 
derivatives  are  connected  to  a,  i,  £2,  and  T  albeit  indirectly. 
Moreover,  there  appears  to  be  no  way  to  elegantly  solve  the 
problem  numerically. 


*_v \v •>  v ’j  ’j  v  vr.vv.v-. 


III.  PARALLAX 


A.  Data  Acquisition 

In  this  instance  the  near-Earth  satellite  was  so  close 
(_<  1000  km)  to  us  that  we  could  separately  resolve  the  two 
streaks  from  the  two  telescopes  at  the  ETS  after  the  videotapes 
were  overlaid  (see  Fig.  3).  The  distance  between  corresponding 
points  on  the  streaks  is  the  parallax.  This  is  related  to  the 
topocentric  distance  of  the  satellite.  Since  we  are  assuming  a 
circular  orbit  and  zenith  passage,  the  topocentric  distance  is 
just  a  -  p  [use  L’Hospital’s  rule  on  the  last  of  Eqs.  (4)  for  an 
analytical  demonstration  of  the  geometrically  obvious]. 

Therefore,  we  can  compute  the  semi-major  axis  a.  The  inclination 
of  the  orbit  comes  from  a  measurement  of  the  slope  of  the 
streak.  Note  that  since  a  is  smaller  for  this  class  of 
satellites  than  those  treated  above,  this  is  an  even  better 
approximation. 

Looking  at  Fig.  3  it  is  imperative  to  realize  that  the 
parallax  is  the  distance  between  corresponding  points  on  the  two 
streaks  and  not  the  perpendicular  distance  between  the  two 
streaks.  We  now  need  to  consider  how  to  go  from  the  experimental 
setup  to  the  appropriate  quantity. 

B.  The  Geometry 

Consider  an  observer  located  at  £  -  p_£  (  x ,<£’).  In  detail, 
this  observer  has  a  geocentric  distance  p,  a  geocentric  latitude 


V.'v. 


v  - 
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PARALLAX 


$ ’ ,  and  a  geocentric  west  longitude  A.  The  local  sidereal  time 
depends  upon  A  and  the  Universal  Time  t  via 


t  -  constant  +  1 .0027379093t  -  A 

The  numerical  constant  exhibited  is  just  t  in  Eq.  (6).  Let  this 
observer  perceive  a  satellite  at  topocentric  position  A  and  A 
when  it  has  geocentric  location  r  -  r£.(a,6).  Equations  (3,  4 , 
and  5)  provide  the  interrelationships. 

Now  consider  a  second,  nearby  observer  whose  geocentric 
location  is  given  by  p  ♦  Ap,  0’  +  A«|>  *  ,  and  A  +  AA.  At  the  same 
instant  of  time  he  will  perceive  the  satellite  to  be  at  A  +  AA 
and  A  +  AA  where 

(p  +  Ap)sin(h  +  Ah) 

tan ( a  -  A  -  AA)  -  - 1 - . - rr- - 7-r-y 

1  -  (p  +  Ap)cos(h  +  Ah) 

(q  ♦  Aq)sin(Y  +  AY  -  6) 

tan(  6  —  A  ~  AA)  “  - 7- - * — ? - -r- - r-rr -r-y 

1  -  (q  +  Aq)cos(Y  +  AY  -  6) 

Although  written  in  a  suggestive  form,  these  equations  are 

exact.  If  |Ap|  and  p  [  ( AAcos<{> ' )  2  +  (A<J>')2]'/2  are  both  small 

quantities  compared  to  p,  then  we  may  simplify  these  results,  for 
AA,  AA,  Ap,  Aq,  and  AY  will  all  be  small  in  magnitude  too. 
Summarizing  several  pages  of  algebra, 


pAXcosh  -  Apsinh  -  p2AX 

A  A  = _ 

1  -  2pcosh  +  p2 

(12) 

-qAYcos(Y  -  5)  -  Aqsin(Y  -  6)  +  q2AY 
1  ~2qco s(Y  -  6 ~+  q* 


AYsecYcscY  =  ( AA/2) { tan[ (a  -  A)/2]  -  tan[h  +  (a  -  A)/2]}  + 
A<J> '  sec<j> '  csc4) '  -  AXtan[h  +  (a  -  A)/2] 


Ap/p  -  Ap/p  -  A4>'tan4>' 

Aq/q  -  Ap/p  +  A^cot^’  -  AYcotY 


and,  since  Ah  =  -AX 

AR/R  =  [ cot ( A  -  Y)  -  cot ( 6  -  Y)]AY  -  AAcot ( A  -  Y) 

C.  Focal  Plane  Formulas 

Once  again  each  observer  records  the  events  as  projected 
onto  his  focal  plane.  If  they  both  have  small  fields  of  view, 
then  the  approximation  in  Eqs.  (9)  is  appropriate.  If  the  first 
observer  measures  5,  n  corresponding  to  A,  A,  then  the  second  one 
will  obtain  £  +  A£  and  n  +  An  corresponding  to  A  +  AA  and  A  +  AA. 
As  AA*  »  -AX  and  AA*  -  A(j 


+  (AA  +  AX ) cosA* 


A£/f  =  -(A  -  A*)A<}>,sinA* 


(13a 


The  separation  between  corresponding  points  on  the  streaks  is 
equal  to  [(A£)2  +  (An)2]1/2. 

However,  all  of  this  occurs  near  a  zenith  passage.  Thus, 
h  =  0,a  =  A  =  A*-T,  6  =  A  =  A  *  -  4> '  ,  Y  -  <j>  ’  ,  and  p  =  q  = 
p/r  -  p/a.  From  Eqs.  (12), 

pAX  -qAY 

A  A  —  - ,  AA— ■ ■»- ; - ,  and  AY — A0  ' 

1  -  p  1  -  q 

Therefore,  Eqs.  (13)  may  be  approximately  written  as 

AXcos<J>'  -A(j>' 

4t/f  •  4n/f  ■  (1 

The  distance  between  corresponding  points  on  the  streak  is 

[  ( AX  cos$  '  ) 2  ♦  (A<|)’)2]1/2 
1  -  p/a 

Were  this  measured  we  could  compute  the  semi-major  axis  a  for  a 
known  observer  separation. 

D.  Data  Analysis  -  Two  Parallel  Lines 

If  all  one  has  is  two  parallel  streaks,  near  the  center  of 
the  field  of  view,  from  two  telescopes  pointing  nearly 


zenithally,  then  one  can  no  longer  determine  which  are 
corresponding  pairs  of  points.  The  reason  is  simple  -  we  are 
looking  at  a  static  picture  wherein  all  time  information 
concerning  the  development  of  the  streaks  has  been  lost.  Hence, 
we  must  proceed  in  another  fashion. 

Let  (£j,ni)  be  the  values  of  5,  n  at  time  t,  for  the  first 
observer.  The  corresponding  values  of  these  coordinates  for  the 
second  observer  will  be  ( £  i  +  +  An).  At  time  t2  the  first 

observer  measures  (£2,n2).  For  the  corresponding  point  on  his 
streak  the  second  observer  will  deduce  (£2  +  A£,n2  +  An).  Note 
that  there  is  no  time  dependence  in  An  or  A£,  see  Eqs.  (14). 

That  is  why  the  corresponding  points  are  identically  offset  and, 
therefore,  why  the  two  streaks  are  parallel. 

Use  the  two-point  form  of  the  equation  of  a  straight  line  to 
deduce  the  equations  for  the  two  streaks.  Then  recast  the  result 
in  the  slope-intercept  form.  The  results  are 

n2  n i  (2^1  -  £  1 n2 

n  »  m£  +  b,  m  -  ■= - =-,  b  -  — ? - = - 

S  2  Si  S  2  SI 

n  -  m5  +  B,  B  -  b  +  An  -  mA£  (15) 

The  two  slopes  are  the  same  but  the  y-intercepts  are  different. 
Clearly  if  we  can  determine  b,  B,  and  m,  then  in  light  of  Eqs. 
(14),  we  can  determine  the  semi-major  axis  a. 


As  our  instruments  and  measuring  devices  are  not  perfect. 


do  not  measure  n  but  approximations  to  them,  say  x,  y. 
Suppose  that  we  make  N  measurements  on  the  fainter  streak 
i(xn,yn)}  and  N  measurements  on  the  brighter  streak 
{(Xn,Yn)l*  Then,  since  it  is  more  difficult  to  measure 
fainter  streaks  than  it  is  to  measure  brighter  ones,  we  weight 
the  measures  on  the  fainter  streak  by  w  e  [ 0 , 1 ] .  A  logical  sum 
of  the  squares  of  the  residuals  to  consider  is 


S  -  I  {w[yn  -  (mxn  +  b)]  +  [Yn  -  (mXn  ♦  B)]  } 

n-1 

Minimizing  S  with  respect  to  m,  b,  and  B  provides  three  normal 
equations  for  their  determination.  From  these  values  and  the 
knowledge  of  p,  A<j> '  ,  and  AX  comes  a  if  the  relationship  (15)  is 
utilized. 

With  a  and  i  determined  and  the  direction  of  motion  known, 
Eqs.  (2)  yields  Mz  and  then  Q  since  we  know  the  time  of 
observation. 


E.  The  Inclination 

There  is  a  purely  geometric  way  to  deduce  the  inclination. 
One  need  only  draw  the  topocentric  celestial  sphere  and  reflect 
upon  the  four  exhaustive  and  exclusive  possibilities  shown  in 
Fig.  4.  By  solving  the  node,  zenith,  and  meridian/equator 
intersection  right  spherical  triangle,  one  comes  to  Eq.  (11b). 


'Ciraai 


NCP 


Figure  4.  Spherical  triangles  for  the  four  exhaustive  and 
mutually  exclusive  kinds  of  zenith  passage. 


Figure  4.  (Continued) 
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